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Abstract 

Fundamental aspects of inverse kinetic theories for the incompressible Navier-Stokes 
equations [Ellero and Tessarotto, 2004, 2005] include the possibility of defining 
uniquely the kinetic equation underlying such models and furthermore, the con- 
struction of a kinetic theory implying also the energy equation. The latter condition 
is consistent with the requirement that fluid fields result classical solutions of the 
fluid equations. These issues appear of potential relevance both from the mathe- 
matical viewpoint and for the physical interpretation of the theory. Purpose of this 
work is to prove that under suitable prescriptions the inverse kinetic theory can be 
determined to satisfy such requirements. 
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1 Introduction 



An aspect of fluid dynamics is represented by the class of so-called inverse 
problems, involving the search of model kinetic theories able to yield a pre- 
scribed complete set of fluid equations advancing in time a suitable set of fluid 
fields. This is achieved by means of suitable velocity-moments of an appropri- 
ate kinetic distribution function /(r, v,t). Among such model theories, special 
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relevance pertains to those in which the state of isothermal incompressible 
fluids is described self-consistently by suitable fluid fields to be considered as 
classical (i.e., strong) solutions of the corresponding fluid equations {regularity 
assumption). In this case the relevant fluid flelds are the mass density, fluid 
velocity and fluid pressure {p, V,p} which are required to be classical solutions 
of the so-called incompressible Navier-Stokes equations (INSE) 



-p+V-(pV) = 0, (1) 

p-^V + Vp + f-//VV = 0, (2) 

V-V = 0, (3) 

pM>0, (4) 

pM>0, (5) 

p(M)=p„>0. (6) 



The flrst three equations (1),(2) and (3), denoting respectively the continuity, 
Navier-Stokes and isochoricity equations, are assumed to be satisfied in the 
open set f2 C (fiuid domain) and in a possibly bounded time interval / C M, 
while the last three inequalities, (4)-(6) apply also in the closure of the fiuid 
domain = ^2 U 5VL. Here the notation is standard. Hence = ^ + V " V 
is the convective derivative, f (r,i) denotes the volume force density acting on 
the fluid element and = ppo > is the constant fluid viscosity, p = [x/ Po 
being the related kinematic viscosity. It is assumed that the fluid flelds and 
f(r,t) are suitably smooth to belong to the functional setting 



VM,p(r,i),f(r,v,t)eCW(nx7), 
< VM,p(r,t)eC(2'i)(Qx/), (7) 
f(r,t) e C(i'°)(n X /), 

[where C^''^\n x /) = C^'^\VL) x C^^)/), with i,j e N]. Hence, {p, V,p} are 
classical solutions of INSE so that the energy equation, obtained by taking 
the scalar product of the Navier-Stokes equation by the fluid velocity V, holds 
identically in the same domain 

^^ + — V- Vp+— V-f-i/V-V^V = 0. (8) 

Dt 2 Po Po 

The set of equations (l)-(6), together with (8), will be denoted in the sequel 
as extended INSE. 

An inverse kinetic theory for INSE is therefore represented by a kinetic equa- 
tion of the form 

i^(F)/ = 0, (9) 
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/(x,t) being a kinetic distribution function, defined in tlie extended pliase 
space r X J, where T = Q x U (with F and U the phase and velocity spaces), 
r = X f/ its closure, while x = (r, v) G F and v G t/ = denote re- 
spectively the state and velocity variables. The distribution /(x,t) is assumed 
suitably regular (i.e., suitably smooth in F x / and summable in F) and to obey 
appropriate initial and boundary conditions. L{F) is an appropriate operator 
(to be denoted as Vlasov streaming operator). This is defined in such a way 
that appropriate velocity moments of (9), are assumed to exist, yield INSE. In 
particular, L(F) can be defined in such a way to allow that the inverse kinetic 
equation (9) admits, as a particular solution, the local Maxwellian distribution 

5/2 

/M(x,t; V,pi) = 3 exp{-Xn. (10) 

(27r)^pf 

Here, the notation is standard [3], thus 



X' = —,, (11) 
< = 2pi/p„ (12) 

Pi being the kinetic pressure. Desirable features of the inverse kinetic the- 
ory involve the requirement that, under suitable assumptions, the functional 
form of the relevant inverse kinetic equation, yielding the INSE equations, be 
uniquely defined. In addition, it might be convenient to impose also the valid- 
ity of additional fluid equations, such for example the energy equation (8). In 
fact, it is well known that the energy equation is not satisfied by weak solutions 
of INSE and, as a consequence, also by certain numerical schemes. Therefore, 
the validity of the inverse kinetic equation yields a necessary condition for the 
existence of classical solutions for INSE. 

Concerning the first point, the prescription of uniqueness on the kinetic equa- 
tion has to be intended in a suitably meaningful sense, i.e., to hold under the 
requirement that the relevant set of fluid equations are fulfilled identically by 
the fluid fields in the extended domain Q x I. This means that arbitrary con- 
tributions in the kinetic equation, which vanish identically under a such an 
hypothesis, can be included in the same kinetic equation. Consistently with 
the previous regularity assumption, here we intend to consider, in particu- 
lar, the requirement that the inverse kinetic equation yields also the energy 
equation (8). 

In a previous work 0, 0|, an explicit solution to INSE has been discovered 
based on a continuous inverse kinetic theory, adopting a Vlasov differential 
kinetic equation defined by a suitable streaming operator L. Basic feature of 
the kinetic equation is that, besides yielding INSE as moment equations, it 
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allows as particular solution local kinetic Maxwellian equilibria for arbitrary 
fluid fluids {po, which belong to the above functional setting (7). How- 
ever, as pointed out in j^, the inverse kinetic equation defined in this way 
results parameter-dependent and hence non-unique, even in the case of local 
Maxwellian kinetic equilibria. This non-uniqueness feature may result as a 
potentially undesirable feature of the mathematical model, since it prevents 
the possible physical interpretation of the theory (in particular, of the mean- 
field force F) and may result inconvenient from the numerical viewpoint since 
|a| may be chosen, for example, arbitrarily large. Hence it would be highly 
desirable to eliminate it from the theory. 

The purpose of this paper is twofold. 

First we intend to prove that under suitable prescriptions the inverse kinetic 
equation can be cast in an unique form, thus eliminating possible parameter- 
dependences in the relevant streaming operator [L(F)].This is achieved by 
analyzing the form of the streaming operator for particular solutions (local 
Maxwellian kinetic equilibria). In this case the the kinetic equation can be 
cast uniquely in an equivalent symmetrized form represented by a Vlasov 
streaming operator. 

As further development of the theory, it is shown that the streaming operator 
can be suitably modified in such a way that the inverse kinetic equation yields 
the extended INSE equations, i.e., besides the incompressible Navier-Stokes 
equations also the energy equation. In particular we intend to prove that the 
mean-field force F can be uniquely defined in such a way that both kinetic 
equilibrium and moment equations yield uniquely such equations. 

The scheme of the presentation is as follows. In Sec. 2 the inverse kinetic equa- 
tion developed in jsf is recalled and the non-uniqueness feature of the mean- 
field force F is analyzed. In Sec. 3 an equivalent representation of the stream- 
ing operator is introduced which permits to define uniquely F. As a result, a 
uniqueness theorem is obtained for the streaming operator L(F). Finally in 
Sec. 4 an extension of the inverse kinetic theory is presented which provides a 
solution also for the energy equation, besides the incompressible Navier-Stokes 
equations. The formulation of the inverse kinetic equation for the extended 
set of fluid equations is obtained by a suitable redefinition of the mean-field 
force F. Also in such a case the vector field F is proven to be unique. 



2 Non-uniqueness of the streaming operator 

Goal of this Section is to investigate the kinetic equation developed in 
to analyze its non-uniqueness features. We start recalling the inverse kinetic 
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equation, which is assumed to be of the form (9) or L{F) f = 0, /(x,t) being the 
normahzed kinetic distribution function associated to the kinetic distribution 
function /(x,t), 

/(x,t) ^ /(x,t)/p,. (13) 
In particular, the streaming operator L is assumed to be reahzed by a differ- 
ential operator of the form 



+ 1; ■ m (15) 

and F(r, v,t; /) an appropriate vector field {mean-field force) defined by Eq.(42) 
(see Appendix A) in terms of vector fields Fq and Fi. As a consequence, 
both Fq and Fiare functionally dependent on the form of the kinetic distribu- 
tion function /(x,t). In particular, requiring that F depends on the minimal 
number of velocity moments (see below), it is defined by Eqs. (42), (43) and 
(44), given in Appendix A. Supplemented with suitable initial and bound- 
ary conditions and subject to suitable smoothness assumptions for the kinetic 
distribution function /(x,t), several important consequences follow j^: 

• the fluid fields {po, V,p} can be identified in the whole fluid domain Q 
with suitable velocity moments (which are assumed to exist) of the kinetic 
distribution function /(x,t) [or equivalent /(x,t)], of the form 

MG{r,t) = Jd\G{^,t)f{^,t), (16) 

where (j(x, t) = l,v,E = ^u'^,vE, uu, and 

u=v - V(r,t) (17) 
is the relative velocity. Thus, we require respectively 



Po = J d^vf{:)c,t), 



V(r,t) = - / rfW(x,t), (19) 

pir,t) = p,ir,t) - Po, (20) 
Pi(r,t) being the scalar kinetic pressure, i.e., 

pi(r,t) = J rfvy/(x,t), (21) 

Requiring, Vp(r,i:) = Vpi(r,t) and Pi(r,t) strictly positive, it follows that 
Po is an arbitrary strictly positive function of time, to be defined so that the 
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physical realizability condition p{r,t) > is satisfied everywhere in Q x I 
(ICR being generally a finite time interval); 

• {po;V,p} are advanced in time by means of the inverse kinetic equation 
Eq.(9); 

• By appropriate choice of the mean-field force F, the moment equations can 
be proven to satisfy identically INSE, and in particular the Poisson equation 
for the fiuid pressure, as well the appropriate initial and boundary conditions 
(see Ref.H); 

• The mean-field force F results, by construction, function only of the ve- 
locity moments (37), i.e., |po, V,pi, Q,2| , to be denoted as extended fluid 
fields .Here Q and n are respectively the relative kinetic energy fiux (de- 
fined in the reference frame locally at rest with respect to the fiuid) and the 
pressure tensor 

Q = J dhn^f, (22) 

U = J d^vuuf; (23) 

• The Maxwellian kinetic distribution function J'm, defined by the equation 
(10), results a particular solution of the inverse kinetic equation (9) if and 
only if {p, V,p} satisfy INSE. 

Let us now prove that the inverse kinetic equation defined above (9) is non- 
unique, even in the particular case of local Maxwellian kinetic equilibria, due 
to the non-uniqueness in the definition of the mean-field force F and the 
streaming operator L{F). In fact, let us introduce the parameter-dependent 
vector field F(a) 



F(a) = F + au ■ VV-aVV ■ u = Fo(a) + Fi (24) 
where F = F(q; = 0), a G M is arbitrary and we have denoted 

Fo(a) = Fo-aAFo = Foa + AiFo(a) , (25) 

AFo = u ■ VV-VV ■ u, (26) 

AiFo(a) = (l + a)u- VV-aVV-u, (27) 



where Fq and Fi given by Eqs.(43),(44). Furthermore, here we have introduced 
also the quantity AiFo(a) to denote the parameter-dependent part of Fo(a). 
In fact, it is immediate to prove the following elementary results: 

a) for arbitrary a G M, the local Maxwellian distribution (10) Jm is a particular 
solution of the inverse kinetic equation (9) if and only if the incompressible 
N-S equations are satisfied; 

b) for arbitrary a in M, the moment equations stemming from the kinetic 
equation (9) coincide with the incompressible N-S equations; 
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c) the parameter a results manifestly functionally independent of the kinetic 
distribution function f{x,t). 

The obvious consequence is that the functional form of the vector field Fq, 
and consequently F, which characterizes the inverse kinetic equation (9) is 
not unique. The non-uniqueness in the contribution Fo(q;) is carried by the 
term aAFo which does not vanish even if the fluid fields are required to satisfy 
identically INSE in the set fl x I. 

We intend to show in the sequel that the value of the parameter a can actually 
be uniquely defined by a suitable prescription on the streaming operator (14) 
and the related mean-field force. 



3 A unique representation 



To resolve the non-uniqueness feature of the functional form of the streaming 
operator (14), due to this parameter dependence, let us now consider again 
the inverse kinetic equation (9). We intend to prove that the mean-field force 
F, and in particular the vector field Fo(a), can be given an unique repre- 
sentation in terms of a suitable set of fiuid fields |po, V,pi, Q,n| defined by 
Eqs. (18)-(21) and (22), (23), by introducing a symmetrization condition on 
the mean field force Fo(q;).To reach this conclusion it is actually sufficient to 
impose that the kinetic energy flux equation results parameter-independent 
and suitably defined. Thus, let us consider the moment equation which corre- 

2 

sponds the kinetic energy flux density G{'x,t) = v^. Requiring that /(x.t) is 
an arbitrary particular solution of the inverse kinetic equation (not necessar- 
ily MaxweUian) for which the corresponding moment q = / d^vv^f (kinetic 
energy flux vector) does not vanish identically, the related moment equation 
takes the form 



-J civG(x,t)/ + V- 1 (ivvG(x,i)/- 



/ dv [Foa + AiFo(a) + F^] • ^^^f 
jd.f 



(28) 



d_ 
di 



G(x,i) + V • VG'(x,t) 



= 0. 



Introducing the velocity moments P2 = J dv^f, £ = / dvuu^f and X — 
J dvuuu/, the kinetic energy flux equation becomes therefore 
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^q + V-P + V- P + 2VQ + VVpi 

at — — 



+ 



+ 



Po 



:Q + Vpi 



D ^ 1 _ 1 
— lnpi + -V- Q-- 

Dt pi Pi 



vn 



•Q - 



2pi 



vn 



1' 


P2_ 


3 


+ 


QV+E 






-2^\ 







2=j 



(29) 



-(1 + a)Q-VV + aVV • Q + -VV : (X + flV 

o — 



2 9V 2 2 

+-— ■ n + -V • VV ■ n+-vv : t 

3 dt = 3 = 3 = 



0. 



Unlike the lower-order moment equations (obtained for (^(x,^) = l,v,'u^/3), 
the kinetic energy flux equation contains contributions which depend linearly 
on the undetermined parameter a. These terms, proportional to the velocity 
gradient VV, yield generally non- vanishing contributions to the rate-of-change 
of q. The sum of all such terms, which are carried respectively by AiFo(a), 
V • VG(x,t) and the convective term V ■ VQ, reads 

M«(/) = 2V- (VQ) - (1 + a)Q-VV + aVV ■ Q. (30) 

which is the contribution to the rate-of-change of q which results proportional 
both to Q (the relative kinetic energy flux) and the velocity gradient (either 
VV or V-V). In order to eliminate the indeterminacy of a, since a cannot 
depend on the kinetic distribution function /, a possible choice is provided by 
the assumption that Mq(/) takes the symmetrized form 

M,(/) = 2V- (VQ) -ivv ■ Q+^Q ■ VV, (31) 

which manifestly implies a = 1/2. Notice that the symmetrization condi- 
tion can also be viewed as a constitutive equation for the rate-of-change of 
the kinetic energy flux vector. In this sense, it is analogous to similar sym- 
metrized constitutive equations adopted in customary approaches to extended 



thermodynamics ll|. On the other hand, Eq.(31) implies Ma(/) = |Qx^, 
^ = V X V being the vorticity field. Thus, Mq,(/) can also be interpreted as 
the rate-of-change of the kinetic energy flux vector Q produced by vorticity 
field ^. From Eq.(31) it follows that Fo(a) reads 



Fo(a = b = --f+ (32) 

^ Po 

+^ (u • VV+VV ■ u) + z/V^V. (33) 
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As a consequence, the functional form of the streaming operator results uniquely 
determined. Finally, for completeness, we notice that the same representa- 
tion for Fo(a) can also be obtained adopting the viewpoint described in 
Appendix B. In fact, since a is functionally independent of the kinetic dis- 
tribution function /(x, t), it can also be defined in such a way to satisfy a 
suitable symmetry condition in velocity-space, which holds in the 'particular 
case /(x, t) = /a/(x, t). This is can be realized by requiring that the Vlasov 
streaming operator L(F) coincides in such a case with a suitable Fokker-Planck 
operator with velocity-independent Fokker-Planck coefficients (see Appendix 
B). 

Finally, it is interesting to point out that the Vlasov streaming operator L{F) 
defined in terms of Fo(q; = |) results by construction Markovian. Hence, the 
position (32) does not confiict with the Pawula theorem 0, 13| which yields 
a sufficient condition for the positivity of the kinetic distribution function /. 
The condition of positivity for the kinetic distribution function satisfying the 
inverse kinetic equation (9) which corresponds to the definition (32) for Fo(a) 
has been investigated elsewhere In particular by assuming that / results 
initially strictly positive and suitably smooth, one can prove that / satisfies an 
H-theorem both for Maxwellian and non-Maxwellian distributions functions. 

As a result of the previous considerations, it is possible to establish the fol- 
lowing uniqueness theorem: 

THEOREM 1 — Uniqueness of the Vlasov streaming operator L{F) 
Let us assume that: 

1) the fluid fields {p, V,p} and volume force density f(r, V,t) belong to the 
functional setting (7); 

2) the operator L{F), defining the inverse kinetic equation (9), has the form 
of the Vlasov streaming operator ( I4 ); 

3) the solution, /(x, t), of the inverse kinetic equation (9) exists, results suit- 
ably smooth in Fx I and its velocity moments |po, V,pi, Q, n| define the fluid 
fields {po,^,p} which are classical solutions of INSE, together with Dirich- 
let boundary conditions and initial conditions. In addition, the inverse kinetic 
equation admits, as particular solution, the local Maxwellian distribution (10); 

4) the mean-field force F{a) is a function only of the extended fluid fields 
|po, V,pi, Q, n| , while the parameter a does not depend functionally on /(x, t); 

5) the vector field AiFo(a) satisfies the the symmetry condition (31). 

Then it follows that the mean-field force F in the inverse kinetic equation (9) 
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is uniquely defined in terms of 



F = Fo + Fi, 



(34) 



where the vector fields Fq and Fi are given by Eqs. (32) and (44)- 



PROOF 



Let us consider first the case in which the distribution function /(x, t) coincides 
with the local Maxwellian distribution /m (10). In this case by definition the 
moments Q, 2 vanish identically while, by construction the mean mean-field 
force is given by F{a) [see Eq.()], a e R being an arbitrary parameter. 

Let us now assume that /(x, t) is non-Maxwellian and that its moment Ma(/) 
defined by Eq.(30) is non- vanishing. In this case the uniqueness of F follows 
from assumptions 4 and 5. In particular the parameter a is uniquely deter- 
mined by the symmetry condition (31) in the moment Ma(/). Since by as- 
sumption a is independent of /(x, t) the result applies to arbitrary distribution 
functions (including the Maxwellian case) . 

Let us now introduce the vector field F'= F+ AF, where the vector field AF 
is assumed to depend functionally on f{x,t) and defined in such a way that: 

A) the kinetic equation L(F')/(x, t) = yields an inverse kinetic theory for 
INSE, satisfying hypotheses 1-5 of the present theorem, and in particular it 
produces the same moment equation of the inverse kinetic equation (9) for 



B) there results identically AF(/m) = 0, i.e., AF vanishes identically in the 
case of a local Maxwellian distribution Jm- 

Let us prove that necessarily AF(/) = also for arbitrary non-Maxwellian 
distributions / which are solutions of the inverse kinetic equation. First we 
notice that from A and B, due to hypotheses 3 and 4, it follows that AF must 
depend linearly on Q, Q — piX.- On the other hand, again due to assumption 
A the vector field AF must give a vanishing contribution to the moments the 
kinetic equation evaluated with respect to G'(x,t) = l,v,E = |-u^. Hence, in 



order that also F' depends only on the moments |po, V,pi, Q,n| (hypothesis 
4) necessarily it must result AF(/) = also for arbitrary non-Maxwellian 
distributions /. 



G'(x,t) = l,Y,E=lu^; 
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4 Fulfillment of the energy equation 



As a further development, let us now impose the additional requirement that 
the inverse kinetic theory yields explicitly also the energy equation (8). 

We intend to show that the kinetic equation fulfilling such a condition can 
be obtained by a unique modification of the mean- field force F = Fo(x,t) + 
Fi(x,i), in particular introducing a suitable new definition of the vector field 
Fi(x,i) [Eq.(44), Appendix A]. The appropriate new representation is found 
to be 



l-V + V- VV + 

at 



F.(x,t;/) = -u|^--V. 

+— f-i/V^V + 

Po 

+— V-Q- ^ fV-nl ■q1+ (35) 
Pi pf L — J J 

As a consequence, the following result holds: 

THEOREM 2 - Inverse kinetic theory for extended INSE 

Let us require that: 

1) assumptions 1-3 of Thm.l are valid; 

2) the mean- field F is defined: 

F = Fo + Fi, (36) 
with Fq and Fi given by Eqs. (32) and (35). 
Then it follows that: 

A) {p, Y,p} are classical solutions of extended INSE in Q x I [equations (1)- 
(6) and (8)] if and only if the Maxwellian distribution function fu (10) is a 
particular solution of the inverse kinetic equation (9); 

B) provided that the solution f{'i<.,t) of the inverse kinetic equation (9) exists in 
r X J and results suitably summable in the velocity space U, so that the moment 
equations of (9) corresponding to the weight-functions G(x,t) = l,v,E = 
exist, they coincide necessarily with extended INSE. 
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C) the two representations (44) ^''^d (35) for Fi coincide identically 



PROOF: 



A) The proof is straightforward. In fact, recalhng Thm.l in [3|, we notice that 
Eqs. (35) and (44) manifestly coincide if and only if the energy equation (8) 
is satisfied identically, i.e., if the fluid fields are solutions of extended INSE. 

B) The first two moment equations corresponding to G'(x,t) = 1, v are man- 
ifestly independent of the form of Fi, both in the case of Maxwellian and 
non-Maxwellian distributions, i.e., (35) and (44). Hence, in such Thm.3 
of applies, i.e., the moment equations yield INSE. Let us consider, in par- 
ticular, the third moment equation corresponding to (^(x,^) = , 



Invoking Eqs. (32) and (35) for Fq and Fi, Eq.(14) reduces to 

piV • V = 

if and only if the energy equation (8) is satisfied. Since by construction pi > 0, 
this yields the isochoricity condition (3). 

C) Finally, since thanks to A) {p, V,p} are necessarily classical solutions of 
INSE, it follows that they fulfill necessarily also the energy equation (8). Hence, 
(44) and (35) coincide identically in Fx/. 

We conclude that (32) and (35) provide a new form of the inverse kinetic 
equation applying also to non-Maxwellian equilibria, which results alternative 
to that given earlier in The new form applies necessarily to classical so- 
lutions. Since weak solutions (and hence possibly also numerical solutions) of 
INSE may not satisfy exactly the energy equation, the present inverse kinetic 
theory based on the new definition given above [see Eq.(44)] for the vector 
field F(x,t) provides a necessary condition for the existence of strong solu- 
tions of INSE. The result seems potentially relevant both from the conceptual 
viewpoint in mathematical research and for numerical applications. 



5 Conclusions 

In this paper the non-uniqueness of the definition of the inverse kinetic equa- 
tion defined by Ellero and Tessarotto (see j^) has been investigated, proving 
that the mean-field force F characterizing such an equation depends on an 
arbitrary real parameter a. To resolve the indeterminacy, a suitably sym- 
metrization condition has been introduced for the kinetic energy flux moment 



d_ 

dt 
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equation. As a consequence, the functional form the mean-field force F which 
characterizes the inverse kinetic equation results uniquely determined. 

Finally, as an additional development, we have shown that, consistently with 
the assumption that the fluid fields are strong solutions of INSE, the mean- 
field force can be expressed in such a way to satisfy explicitly also the energy 
equation. 

The result appears significant from the mathematical viewpoint, the physical 
interpretation of the theory and potential applications to the investigation of 
complex fluids, such as for example those treated in 0,13). In fact, it proves 
that the inverse kinetic theory developed in can be given an unique form 
which apphes to classical solutions of INSE. 
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6 Appendix A: relevant velocity moments and mean-field force 



Here we recall for completeness the expressions of the velocity-moments of the 
kinetic distribution function and of the mean-field force F given in The 
relevant moments are of the form 

M{G) = J dhG{^,t)f{^,t), (38) 

G'(x, t) being the weight functions G(x, t) = 1, v,m^/3,u u^/3, u u. In partic- 
ular, we identify the fluid fields {p = po, V,p} with 

Po = J dhf{^,t), (39) 

V(r,t) = - / ciW(x,t), (40) 

p{r,t)=pi{r,t)-Po, (41) 

Po being a strictly positive real constant defined so that the physical realiz- 
ability condition p{r,t) > is satisfied everywhere in the closure of the fluid 
domain Q. Finally the expression of the mean-field force reads: 

F(x,t;/)=Fo(x,t;/) + Fi(x,t;/), (42) 
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where Fq and Fi are the vector fields: 



Fo(x,t;/) 



— v-n-vpi-f +u-w+uv^y, 

Po^ — ^ 



(43) 




Q= 0. 



(45) 
(46) 



7 Appendix B: Fokker-Planck representation 

It is interesting to point out that the choice of the parameter a = 1/2 and 
of the mean-field force Fo(q;) (32) can be obtained also by requiring that the 
Vlasov streaming operator (14) results "equivalent" to an appropriate Fokker- 
Planck operator (i.e., that it yields an inverse kinetic theory for INSE which 
admits local Maxwellian equilibria). 

Since by assumption the parameter a does not depend functionally on /(x, t) it 
is sufficient to impose its validity only in the case of local Maxwellian equilibria 
/m- In such a case, it is possible to require that there results identically 



^F)/m= Ut + v~ /m+ E -77r-iWM} = LMF)fM, (47) 



where the Fokker-Planck coefficients Dj are assumed velocity-independent and 
Lfp{F) denotes a Fokker-Planck operator which by construction is equivalent 
to the Vlasov operator L{F). There results 





Di(r,i;/M) 



-— f + i/V'V, 

Po 



(48) 




(49) 



14 



M, 



(50) 



On the other hand, the relationship (47) holds if and only if 



F(«) = Di-^ uR^+U^ u 



Pi 



1 uu 

2^ + 4— 



■ D (r,t), (51) 



namely the parameter a necessarily results equal to 1/2, i.e., 



Fo(« = I) 



— f + ^ (u ■ VV+VV ■ u) + z/V^V. 
Po 2 



(52) 
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